Berry phases of quantum trajectories in semiconductors under strong terahertz fields 
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Quantum evolution of particles under strong fields can be essentially captured by a small number of quantum 
trajectories that satisfy the stationary phase condition in the Dirac-Feynmann path integrals. The quantum trajec- 
tories are the key concept to understand extreme nonlinear optical phenomena, such as high-order harmonic gen- 
eration (HHG), above-threshold ionization (ATI), and high-order terahertz sideband generation (HSG). While 
HHG and ATI have been mostly studied in atoms and molecules, the HSG in semiconductors can have interest- 
ing effects due to possible nontrivial "vacuum" states of band materials. We find that in a semiconductor with 
non-vanishing Berry curvature in its energy bands, the cyclic quantum trajectories of an electron-hole pair under 
a strong terahertz field can accumulate Berry phases. Taking monolayer M0S2 as a model system, we show that 
the Berry phases appear as the Faraday rotation angles of the pulse emission from the material under short-pulse 
excitation. This finding reveals an interesting transport effect in the extreme nonlinear optics regime. 
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After excitation by a weak near-resonant laser in semicon- 
ductors, an electron and a hole can be created and driven 
into large amplitude oscillations by an intense low-frequency 
ac electric field such as from a terahertz (THz) laser. The 
recollisions between the electron and hole during the os- 
cillations will generate high-order sidebands relative to the 
excitation frequency |T||2). This high-order THz sideband 
generation (HSG) has potential electro-optical applications 
such as wide-band optical multiplexers, terabit/sec optical 
communications, and optical pulses with ultra-high repetition 
rate EIH. The HSG is analogous to the high-order harmonic 
generation (HHG) in atoms and molecules BUS). Both HSG 
and HHG spectra are characterized by a wide-band plateau 
with a sharp cut-off. The HSG spectrum has been well under- 
stood using the quantum trajectory theory developed in HHG, 
in which the quantum evolution of particles under a strong 
field is described by a few paths that satisfy the stationary 
phase condition (i.e., the saddle points) in the formalism of 
Dirac-Feynmann path integrals HUH. The HSG or HHG cut- 
off is determined by the maximum energy the particles in the 
quantum trajectories can acquire from the driving field HIE]. 

A fundamental difference between HHG in atomic systems 
and HSG in semiconductors is that the "vacuum" state of 
a semiconductor can have non-trivial structures (such as in 
topological insulators EG)). When the eletron (or hole) in a 
spin-orbit-coupled semiconductor is accelerated by an ac elec- 
tric field F (f), not only does the quasi-momentum evolve ac- 
cording to the semiclassical equation k = -eF (f) [8 1, but also 
its Bloch wavefunction (the direction of the spin) is changed. 
Thus the evolution driven by the electric field leads to a ge- 
ometric phase in addition to the dynamical one, which is the 
famous Berry phase in a cyclic evolution J9). This geomet- 
ric phase is of fundamental importance for a gauge-invariant 
description of the nonlinear optics in insulators IfTOifTSl . 

The Berry phase effect is clearly seen in the representation 
of the polarization operator in the basis of Bloch states, which 



is found by Blount Ifl3l 

(</f„,k| r \lf/ m , k >) = i [dnmVk + (U n ,k\ V k |«m,k)] (5 (k - k') , (1) 

where the first term gives change of the quasi-momentum and 
the second term is the so-called Berry connection or Berry 
vector potential. Through the polarization operator, the Berry 
phase (or, more intrinsically, the Berry curvature) appears nat- 
urally in various optical effects in condensed matter systems 
as revealed by some recent works. For example, the presence 
of the Berry phase effect was noticed in the optical birefrin- 
gence effects of a pure spin current [ 14 1 or the second-order 
non-linear spectroscopy of spin currents ifTTl . Also the photo- 
galvanic effect in topological insulator surfaces can depend on 
the Berry curvature |[T5l . The interaction between an intense 
THz laser and semiconductors |2] provides a new opportu- 
nity to explore the Berry phase effect in the regime of extreme 
nonlinear optics. 

In this Letter, we show that the optical response of a semi- 
conductor under an intense THz field explicitly includes the 
Berry phase. We analyze the effect using the quantum tra- 
jectory theory and apply the theory to monolayer M0S2 as a 
model system. In the time-domain response, we find that the 
Faraday rotation angle of the emission delayed by integer mul- 
tiples of the THz laser period is given by the Berry phase of 
a specific trajectory. The quantum trajectory approximation is 
verified by numerical simulations. 

Let us consider a general semiconductor under a strong 
THz field F(f), which enters into the Hamiltonian through a 
uniform electromagnetism vector potential: p — > p + eA (f), 
with F = -dA/dt. Because A preserves the translational 
symmetry, Bloch's theorem still applies and we write the 
Hamiltonian in the k-space representation ff (£(?)) fl6l . with 

k (f) = k + eA (f). The instantaneous Bloch states of H (k (f)) 
are obtained from the original Bloch states by simply chang- 
ing k to k, 

H(k(t))\±,»Mi)) = E±J±,»Mi)), (2) 
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where + and - are the indices of the conduction and valence 
bands, respectively, and fi is the spin index introduced to indi- 
cate possible band degeneracy of the system. 

Now let us calculate the linear response of this system to 
a near infrared (NIR) laser that creates electron-hole pairs 
at the band edge with interaction Hamiltonian 77nir = _ P ■ 
ENiRe~' n ' + h.c, where the interband polarization operator P 
in the interaction picture is 



P(0 



dke 1 .ti 



- k ^V,k (0 • 



(3) 



Here e and h are electron and hole operators, respectively, and 
^jv,k (0 = -ie (<A+,/j,k (0| |iA-,v,k (0) i s tne interband dipole 
moment, with |iA±,^,k (0) denoting the adiabatic evolution of 
the instantaneous Bloch states under the driving of the THz 
field 



\tfr±jiM W) = |±.«,k(0) 



Te 



-i f Bi dr+i f 

J-co k(r) J-a 



(4) 



where t is the time-ordering operator, the Berry connection 
is defined as if^^) = * (±,/*,k| V k |±, v,ky, and summation 
of repeated dummy indices is assumed. In general, the Berry 
connection can be non-Abelian. Assuming the initial state is 
the vacuum state |G) with empty conduction bands and filled 
valence bands, we obtain the linear response to the NIR opti- 
cal field as 



<P(0> 



= i I dt' J — 

J-co J (2n 



^fi'v'Mn ■ e nir 



dt' (G\P(t)H mR (t')\G) 



Te 



(5) 



where d is the dimension of the system, e k 

is the energy of the electron-hole pair and d 

-ie(+, y,k| V k |— is the instantaneous dipole moment. 
The Berry phase enters Eq. |5]) due to the requirement of the 
gauge invariance of the physical result under the local gauge 
transformation |±,/i, k) — > |±,ar,k) (k), which introduces 
the gauge freedom of d VjU k . For the case without band degen- 
eracy, the Berry phase becomes Abelian, and the response is 
reduced to 

<P(0> = (irf I-J f f ' Emr 



T)' 



(6) 



where ^4 = -^C 



is the combined Berry connection of 



the electron-hole pair. It is interesting to note that the Berry 
term ^ ■ k in the phase factor can be regarded as a general- 
ization of the interaction energy A • r in the electromagnetism 
theory in real space to the momentum space fl6l . 

Now let us focus on the geometric phase part of Eq. (|6j. 
Without loss of generality, we consider an elliptically polar- 
ized THz field in the x-y plane 



F (0 = F (cos (6) cos (cjt) , sin (6») sin (cat) , 0) . 



(7) 



Then the electron-hole pair goes along an elliptical path in the 
k-space under the driving of the THz field: 

k (f) = (k x - Icq cos sin (a>t) , k y + &o sin cos (cot), , (8) 

where Icq = eF/a>. After an integer multiple periods of 
the THz field t — f = nT = Inn/a), the electron-hole pair 
completes a closed path in the k-space. During this cyclic 
evolution, the geometric phase acquired by the electron-hole 
pair equals the Berry curvature times the area S (9, ko) = 
nnkZ sin (26) 1 2 enclosed by the path, which is further related 
to the polarization and strength of the THz light. 

To measure the Berry phase effect of quantum trajectories, 
we can apply a short NIR laser pulse to the system at time 
t' — to, with the width of the pulse much smaller than T. Thus 
the pulse can be approximated by a 5-pulse Enir = E<5 (t — to). 
Then response of the system at t„ = to + nT explicitly contains 
the Berry phase of a closed path: (p£' (k) = J" stf^ ■ iik(r). 
In order to separate the Berry phase from the dynamical phase, 
we can introduce some specific interference that singles out 
the geometric phase part. We note that under the time-reversal 
transformation, the direction of the path in k-space is reversed 
and the Berry phase becomes opposite, while the dynamical 
phase is unchanged. It leads us to consider the solid state sys- 
tems that preserve the time-reversal and inversion symmetry 
and have nontrivial Berry phases, such as the topological in- 
sulators ||6]|7), monolayer M0S2 and other group- VI dichalco- 
genides [17-21] and bilayer graphene l22l l23l . We denote 
one state of the Kramers pair as the pseudospin state ft and the 
other as JJ.. From the time-reversal and inversion symmetry, we 
obtain the following relations k = dJn_ k = _< ^J| Jik := d k , 
%k = %,-k := Sk and 04)^ = (^k)^j := ^4- Thus we get 
the key formula 



^(k) 



(k) 



(k) 



f 

J to 



stt ■ dk, 



(9) 



i.e. the Berry phase of the two time-reversal related paths 
are opposite to each other. With these considerations, the re- 
sponse at t n is simplified as 



<P 



+i f ^ e " ! ' xr " (T) " 4 " )(k)d ^) d k)- E - ^ 



Here (P(f„)) is given by the interference between two kinds 
of responses with opposite Berry phases. 



Equation ( 10 1 can be studied using the quantum trajectory 
theory, i.e. the stationary phase formalism flT] [5] . In the path 
integral, the electron-hole pairs move along all possible tra- 
jectories when driven by the THz field, with the phase given 
by the action S ± (k) = JT " (e k(T) ± ■ eF (t)) dr. As the THz 
field is strong, the motion amplitude of the electron-hole pair 
is much larger than the quantum fluctuation. Thus the re- 
sponse is dominated by the stationary phase points of the ac- 
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tions 



V k S ± (k)= f l (V k s i(T) ±V k ^ i -e¥(T))dT = 0, (11) 

J to 

plus the Gaussian quantum fluctuation around them. V k e k - 
v k is the semiclassical velocity of the electron-hole pair. 
The Berry connection term gives a gauge dependent motion 
V k (affc ■ e¥) = V k • e¥) - (e¥ ■ V k ) ^ - Be- 
cause we are considering the cyclic evolution along a closed 
loop, the last gauge dependent term vanishes and the first two 
terms gives a gauge-invariant physical quantity Q^eFye,-, 



with Q 



k.k. 



dk being the Berry curvature. We 



note that this is the well-known anomalous velocity that is 
responsible for various Hall effects lfl6l l24l |251 . The station- 
ary phase condition in Eq. (Hi therefore means the return of 



the electron to the hole after nT under the acceleration by the 
THz field. The Berry phase (pjj (k) in the actions is gener- 
ally a slowly-varying function of k and much smaller than the 
dynamical phase factor. Therefore the stationary phase points 
are determined by J f " v^dr = 0, which has a simple solution 
v k = if the effective mass model is used. This means that the 
response is dominated by the trajectories of the electron-hole 
pairs whose paths in the k-space are centered at the extreme 
points of the energy band (see Fig. [TJ. Thus Eq. (jTOj is ap- 
proximated by 



<PfeO> 



2ldk P ~' C s feT) dT <R \a* a . . F l 



To be specific, from now on we consider the monolayer 
M0S2 as the model system. This material has interesting spin- 
valley coupling and has potential applications for novel spin- 
and valley-based information processing lfT7l42TI . The effec- 
tive Hamiltonian describing the Bloch states at the band edges 
is given by lfT8ll 

H (k) = A (t1c x o- x + k y o- y ) + Mcr z , (13) 

where 2M = 1.59 eV is the band gap, A = 3.51 eV • A 
and t = ±1 is the index of the +K valley. The energy spec- 
trum is e k = 2 VM 2 + A 2 k 2 and the stationary phase point is 
k = 0. The two valleys are related by time-reversal trans- 
formation and we denote the state at the —K valley as the 
pseudospin state f|\ We choose the gauge such that the Berry 
connection ^4 = E ^, M (k x e y - k y e x f and the dipole moment 

d k ~ d cv x (e x + !'e y ) for small k in the —K valley. For the 
sake of simplicity we have not included the Coulomb inter- 
action between electrons and holes. This is justified for the 
band edge excitation since the exciton binding energy (100s 
of meV ||26l ) is much greater than the THz field and therefore 
the exciton bound states are far off-resonant from the NIR ex- 
citation. 
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FIG. 1 : (color online). Schematics of optical response to a short NIR 
pulse and quantum trajectories in semiconductors under strong THz 
fields, (a) An electron-hole pair is excited by an NIR pulse and then 
driven along quantum trajectories by an elliptically polarized THz 
field. At t„ — tg = nT, the electron-hole pair completes an elliptical 
loop in k-space (see (b)) and obtains a nonzero Berry phase. The 
polarization of optical emission at t„ is rotated by an angle nip that 
is equal to the Berry phase ipg . (b) Paths of the electron-hole pairs 
in k-space, where the green solid curve is for the electron-hole pair 
satisfying the stationary phase condition in Eq. i | 1 1 fr . The red dashed 
one does not satisfy Eq. l[TTJ. (c) The green solid curve gives the 
quantum trajectory of an electron-hole pair that recombines at t\ . The 
red dashed curve corresponds to the red dashed path in (b), in which 
the electron-hole pair does not make a close path in real space. 



We assume that the NIR field is linearly polarized in the x-y 
plane with E = Ee\\. The dipole moment gives 

d k d k = \d cv , k | 2 [(e. v e t + e T e v ) + i (e x e y - e y e x )] . (14) 

The real part of ( [T4| leads to the longitudinal response along ey 
while the imaginary part leads to the transverse response along 
e_L (which is related to the Faraday rotation of the emission) 
with ey x e ± = e z . Thus the longitudinal response of Eq. ( 12 1 
is 

<P( f „))„ =cos«>) k=0 j ^^^KmJe, 



(15) 



and the transverse response is 



(16) 

The Faraday rotation angle of the optical emission is exactly 
given by the Berry phase 



0« (0,fc Q ) = ^(k = O)*^° 



sin (20), (17) 
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where A 2 /M 2 is the Berry curvature for small k ifTSll . 

The Faraday rotation effect of an elliptically polarized THz 
field can be intuitively understood as illustrated in Fig. |2ja). 
The electron-hole pair created by a linearly polarized short 
NIR pulse excitation is a superposition of the valley states 
|ft) + III). After the cyclic evolution under the THz field, the 
states at +K valleys (corresponding to r = ±1) obtain the 
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FIG. 2: (color online). Faraday rotation of the optical emission of 
monolayer M0S2 under a strong THz field, (a) Schematic for op- 
posite Berry phases of quantum trajectories in different valleys, (b) 
and (c) plot Faraday rotation angle of the optical emission at 
t = nT (n = 1,2, 3) as a function of the polarization ellipticity (6) 
and the strength (F) of the THz field, respectively. The symbols 
show the numerical integration (Nl) results and the lines are the sta- 
tionary phase approximation (SPA) in Eq. \17\ . In (b), the THz field 
strength F = 8 kV/cm (i.e., k = 0.02). In (c), the THz field is 
circularly polarized (i.e., 9 = n/4). 



same dynamical phase (po and opposite Berry phases +<ps, 
which are the Berry curvature fluxes through the area en- 
closed by the quantum trajectories. Thus the final state is 
e "t>D ^ g i0s ^ + e -ups ijj^ w hich results in emission with linear 
polarization rotated by an angle <p B . 



Equation ( 17 1 shows clearly the effect of the Berry phase 
on the optical response of the semiconductor in an intense 
THz field and provides a new method to directly measure the 
Berry phase of the energy bands in momentum space. Since 
the Faraday rotation angle is independent of the strength of 
the optical response, the measurement does not rely on the 
specific form of the energy spectrum, the value of c/ c ,,,k, or 
the dephasing of the electron-hole pair during the evolution. 
However, the result does depend on the optical selection rule 
of the dipole moment dk ~ (e A + /e v ), which is due to the rota- 
tional symmetry of the system. Thus Eq. ( [T7} may be applied 
to other two-dimensional spin-orbit coupled semiconductors 
with (approximate) rotational symmetry. For a semiconductor 
having a different form of the dipole moment, the Faraday ro- 
tation angle of the optical emission may have a more complex 
relation to the Berry phase. 

To verify the validity of Eq. ( fT7| ), we compare it with the 
numerical results obtained from the standard numerical inte- 
gration of Eq. ([6]) for a M0S2 monolayer. The frequency of the 
THz field is u> = 4 meV and the NIR pulse has the gaussian 
form £e A exp where Q = 2M and the width of 

the pulse is such that todt = 0.2 («: 2n). Some results calcu- 
lated for different n, 9 and field strength F (i.e. £0) are shown 
in Fig. [2] where the lines are the stationary phase approxima- 
tion and the symbols give the numerical integration results. 
We can see that Eq. ( 17 1 is a good approximation. 

In the discussions above, we only consider the emissions 
at t — fo + nT, However, in the case of t - to + nT, the 
trajectories of the electron-hole pairs also obtain a geometric 
phase. For the monolayer M0S2, the dipole moment is nearly 
constant for small k in the gauge we chose above. Then based 
on the same reasoning, we see that the Faraday rotation angle 
of the emission at any time t is given by the geometric phase of 
the trajectory that satisfies the stationary phase point equation 
£vk (T A = 0. 

In summary, we have obtained a Berry phase dependent 
theory of optical response in spin-orbit-coupled semiconduc- 
tors under strong THz fields, where the Berry phase enters 
the formula as required by the gauge invariance of the opti- 
cal response. This theory is investigated using the quantum 
trajectory theory and applied to the monolayer MoS2- The 
Faraday rotation angle of the optical emission is exactly equal 
to the Berry phase of the quantum trajectory that satisfies the 
stationary phase condition. This result can be generalized to 
semiconductors without time-reversal symmetry. Even more 
interesting, the theory can be applied to semiconductors with 
non-Abelian Berry connection such as the three-dimensional 
topological insulators |6||7). The quantum trajectory will then 
be accompanied by a nontrivial (pseudo)spin rotation that is 
determined by the non-Abelian Berry phase. 



5 



This work is supported by Hong Kong RGC/GRF 401512 
and the CUHK Focused Investments Scheme. 



* Electronic address: rbliu@phy.cuhk.edu.hk 

[1] R.-B. Liu and B.-F. Zhu, in Physics of Semiconductors: 28th 
International Conference on the Physics of Semiconductors - 
ICPS 2006, edited by W. Jantsch and F. Schaffler (AIP, New 
York, 2007), vol. 893 of AIP Conf. Proc, pp. 1455-1456. 

[2] B. Zaks, R. B. Liu, and M. S. Sherwin, Nature 483, 580 (2012). 

[3] J. L. Krause, K. J. Schafer, and K. C. Kulander, Phys. Rev. Lett. 
68, 3535 (1992). 

[4] P. B. Corkum, Phys. Rev. Lett. 71, 1994 (1993). 

[5] M. Lewenstein, P. Balcou, M. Y. Ivanov, A. LHuillier, and P. B. 
Corkum, Phys. Rev. A 49, 2117 (1994). 

[6] X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 (201 1). 

[7] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010). 

[8] N. W. Ashcroft and N. D. Mermin, Solid State Physics (Saun- 
ders, Philadelphia, 1976). 

[9] M. V. Berry, Proc. R. Soc. Lond., Ser. A 392, 45 (1984). 
[10] J. E. Sipe and A. I. Shkrebtii, Phys. Rev. B 61, 5337 (2000). 
[11] J. Wang, B.-F. Zhu, and R.-B. Liu, Phys. Rev. Lett. 104, 256601 
(2010). 



[12] K. S. Virk and J. E. Sipe, Phys. Rev. Lett. 107, 120403 (2011). 
[13] E. I. Blount, in Advances in Research and Applications, edited 

by F. Seitz and D. Turnbull (Academic Press, New York, 1962), 

vol. 13 of Solid State Physics, pp. 305-373. 
[14] J. Wang, B.-F. Zhu, and R.-B. Liu, Phys. Rev. Lett. 100, 086603 

(2008). 

[15] P. Hosur, Phys. Rev. B 83, 035309 (201 1). 
[16] D. Xiao, M.-C. Chang, and Q. Niu, Rev. Mod. Phys. 82, 1959 
(2010). 

[17] K. F. Mak, C. Lee, J. Hone, J. Shan, and T. F. Heinz, Phys. Rev. 

Lett. 105, 136805 (2010). 
[18] D. Xiao, G.-B. Liu, W. Feng, X. Xu, and W. Yao, Phys. Rev. 

Lett. 108, 196802 (2012). 
[19] H. Zeng, J. Dai, W. Yao, D. Xiao, and X. Cui, Nature Nanotech. 

7, 490 (2012). 

[20] K. F. Mak, K. He, J. Shan, and T. F. Heinz, Nature Nanotech. 7, 
494 (2012). 

[21] T. Cao, G. Wang, W. Han, H. Ye, C. Zhu, J. Shi, Q. Niu, P. Tan, 
E. Wang, B. Liu, et al., Nat Commun. 3, 887 (2012). 

[22] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, 
and A. K. Geim, Rev. Mod. Phys. 81, 109 (2009). 

[23] E. McCann and V. I. Fal'ko, Phys. Rev. Lett. 96, 086805 (2006). 

[24] M.-C. Chang and Q. Niu, Phys. Rev. Lett. 75, 1348 (1995). 

[25] G. Sundaram and Q. Niu, Phys. Rev. B 59, 14915 (1999). 

[26] X. Xu, private communication. 



